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Abstract 

Influence of volume capture on a process of volume reflection of ultrarelativistic par- 
ticles moving in bent single crystals was considered analytically. Relations describing 
various distributions of particles involving in the process, the probability of volume cap- 
ture and behavior of channeling and dechanneling fractions of a beam were obtained. 
Results of study will be useful at creation of multicrystal devices for collimation and 
extraction of beams on modern and future accelerators. 



1 Introduction 

Volume reflection of charged particles in single crystals represents the coherent scattering 
of these particles by planar or axial electric fields of bent crystallographic structures. For 
the first time, this effect was predicted in Ref. [1] on the basis of Monte Carlo calculations. 
Recently in Ref. [2] an analytical description of volume reflection of ultrarelativistic par- 
ticles was considered. Besides, this process was observed and investigated in a number of 
experiments [3, 4, 5]. Results of study of volume reflection is interest to series of important 
applications for accelerator technics (as beam collimation, extraction and others). For work 
of these applications is assumed to use many times repeated volume reflection process ( in 
special crystal structures). One can expect that due to a high efficiency of the process the 
sum efficiency of crystal devices will be also high enough. 

The main process which influents on an efficiency of volume reflection is volume capture 
[6] when moving in bent single crystals particles may be (due to multiple scattering on 
atoms) captured in a channeling regime. The analytical description of volume reflection (see 
[2]) gives good predictions for many characteristic of the process, but volume capture was 
not considered in this theory. 

This paper is devoted to complete analytical consideration of volume reflection, including 
influence of volume capture on the process. Main results of Ref. [2] remain valid but their 
simple modernization is required. No doubt that Monte Carlo simulations give detailed infor- 
mation about process at some specific initial conditions. However, an analytical description 
allows one to make observation of the problem as a whole. 

Our experience [7] shows that an accuracy of the analytical calculations of such param- 
eters as a mean and mean square angles at volume reflection compares with the results of 
Monte Carlo simulations. Besides, in Ref. [2] is shown how to use generalized parameters for 
description of the process. Because of this, one can obtain set of results from one calculated 
result. It is important for emulation of the process, when the experimental results are known 
at some conditions (for example, for one energy of beam) and these results should be used 
at another conditions (another energy). 

The paper is organized as follows. First, we describe our method for study of diffusion 
processes in bent single crystals. Than based on this method we find the probability of 
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volume capture and energy distribution underbarrier particles. After that we investigate 
behavior of the channeling fraction at its propagation in single crystal. As result we get the 
analytical relations for distribution functions of particles. 

In the second part we illustrate the obtained relations by numerical calculations and then 
we discuss and summarize presented in the paper results. 



2 One-trajectory approximation of a diffusion process 

Let us consider motion of the charged ultrarelativistic particle in bent single crystal. In the 
absence of multiple scattering we can write the following equation [2]: 

E = E [3 2 x 2 /(2c 2 ) + U(x) + (3 2 E x/R, (1) 

where E is the particle energy, E is the transversal energy, U(x) is the periodic interplanar 
potential as a function of the transversal coordinate x, R is the bending radius of a single 
crystal, (3 is the particle velocity divided on the velocity of light c. 

Obviously, this equation reflects the conservation of the transversal energy. Introduction 
of multiple scattering violates this condition. Let us suppose that the moving in bent crystal 
one separate particle is scattered on the angle A9 (by multiple scattering process). Then 
one can get 

AE = E (3 2 9 m A9 + ^E (A9) 2 . (2) 

Here m ~ x/c is the angle of particle due to regular motion in accordance with Eq.(l) 
Averaging over different values of AO we obtain the following equations 

AE = \e {A6) 2 , (3) 

(E' — E — E) 2 = E 2 (3 2 6 2 m (A6) 2 = A(E - U(x) - (3 2 E x/R)AE (4) 

where E' , E are the variable value and mean energy losses of the transversal energy. 

After averaging over different angles for different particles we get distribution function 
over changing transversal energy (from distribution function over plane scattered angle): 

F(E', I) = -=L- exp -[(Ef -E- E) 2 /(2a 2 E )}, (5) 



where E and a E are expressed through the mean square angle of multiple scattering 9 2 : 

E = l -E Q 6\ (6) 
a E = E 2 [3 2 9 2 Ji. (7) 
Note that function F satisfies to the following equation: 

OF _dE OF lda 2 E d 2 F 

Thus we got simple distribution function for one trajectory. This function depends on only 
two parameters: E and a E . These parameters are independent of the current transversal 
energy E . It is clear that obtained here approximation is valid on a short distance, when 
parameters E and a E are small enough with comparing of the characteristic value SE. 

More traditional consideration of a similar problem is based on the diffusion equation with 
the diffusion coefficients depending on the transversal energy. This fact strongly complicate 
possibility for solution. 
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3 Volume capture probability 



Fig. 1 illustrates behavior of the effective potential U e (x) — U(x) + /3 2 Eqx/R as a function 
of the transversal coordinate x. The main problem of our consideration is determination of 
an inefficiency of volume reflection of charged particles in bent single crystals. Mechanism of 
the inefficiency is connected with the volume capture process in these structures. Really, the 
moving particle has some probability to lose (because of multiple scattering) its transversal 
energy and to continue motion in the channeling regime. However, the captured particle 
has some probability (due to dechanneling) quickly return in underbarrier motion and to be 
similar as usual particle at volume reflection. This fact requires the detailed description of 
the behavior of volume captured particles. 

In this section we find the probability of volume capture. Fig. 1 illustrates this pro- 
cess. Obviously that volume capture takes place mainly in the vicinity of critical (reterned) 
point. In Fig. 1 this region corresponds to x\ < x < x c for transversal energies of par- 
ticle E m < E < E M . In Ref.[2] was shown that particles of beam with the small angle 
divergence in absence of multiple scattering are equiprobably distributed over transversal 
energies. Multiple scattering changes the initial distribution of particle over transversal en- 
ergies. However, the particle distribution in every energy period SE keeps approximately 
equiprobable. Evidence is similar to considered case in Ref. [2]. 

Multiple scattering we consider as in an amorphous medium but instead constant nuclear 
density we use variable one in accordance with a current particle position. For this one can 
give different theoretical likelihood explanations but more evident argument is fact that 
Monte Carlo calculations at the above mentioned assumptions give a good agreement with 
experiments. 

Taking into account stated assumptions we define the probability of volume capture as an 
normalized number of particles (with the initial energy range from E m till E M = E m + SE) 
which obtain (due to multiple scattering) the transversal energy corresponding to finite 
planar motion (channeling). 

As illustrated in Fig. 1 for every particle with the transversal energy E m < E < Em 
and coming to coordinate x\ there are three different possibilities (see numbers near curves 
in Fig. I): l)to undergo reflection; 2)to lose the transversal energy and to become particle 
captured in channeling regime; 3)to increase the transversal energy, so that E > E M . In 
the case of number 3 theoretically the similar three possibilities is conserved but for small 
enough bending radii the third variant is impossible practically. At condition of smallness 
of bending radius the probability for particle (with E m < E < Em) to lose significantly 
the transversal energy before coordinate x 1 is very small. Thus, we see that the volume 
capture probability one can represents as a sum of the two terms E\ = e ia + 6^, where Si a 
is the probability of volume capture in the coordinate range from x\ till X2 and e\b is the 
probability in the coordinate range from x<i till x^. 

Let us calculate the first term S\ a . Our consideration will be based on one-trajectory ap- 
proximation. Eq.(5) gives the distribution for one trajectory of particle with the transversal 
energy E. The summary distribution function is a result of integration over all the range of 
transversal energies (from E m till Em)- Taking this into account we get the probability to 
change the transversal energy E of a particle less than the critical energy E c in the following 
form: 



where £ m (£) = E/5E and £(£) = cr E /5E, E,a E are the mean and square mean losses of 
transversal energy (see Eqs.(6)-(7)), 5E = E /3 2 d/R, erf(x) = 2/0F f£ exp(-t 2 )dt. One can 




(9) 
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write the following equations for the mean and mean square energy losses: 



E = 2A 2 (L /X )/E 2 , (10) 



a E = A^/2(L 1 /X )/E . (11) 

Here A ?a 10 — UMeV is the constant value, X is the radiation length of a single crystal 
and L and L\ functions are: 

L = f°£ f\ Pa {x(t)) + p e {x(t))/Z 2 )dt ) (12) 
2p Jti 



Ll = f t2 (p a (x(t)) + Pe (x(t))/Z 2 )(E - U(x(t)) - E (3 2 x(t)/R)dt, (13) 



2E c 
Po -Ai 

where p a (x),p e (x) are the planar atomic and electron densities for a selected plane of single 
crystal, Z is the atomic number, p ~ Nq/Vq is the mean atomic density (iVo is number atoms 
in a fundamental cell of volume Vq). t\ is time which corresponds to finding of particle in 
the nearest local maximum of U(x) before critical point x c (for this particle). Time t 2 
corresponds to location of the particle in a critical point. The coefficient 2 in Eqs. (11) and 
(12) is due to symmetry of particle motion before and after a critical point. 
Taking into account the relation: 

dt = . , X =, (14) 

^[E-U(x)-E (3 2 x/R] 

we can represent Eqs. (11) and (12) in the following form: 

e = * r ±£M±mm dx , m 

V2Eop X Jxi ^ E -U(x) - E Q f3 2 x/R 



2i A 

= 1 



[ r (p a (x) + p e (x)/Z 2 )^E - U(x) - E (3 2 x/Rdx]l (16) 



E£(p X )2 

It is convenient to rewrite Eqs. (15) and (16) in the following more general form: 
A 2 d rvc (p a (y)+p e (y)/Z 2 ) , A 2 d 



E 



_ . r (PM+mm iy . Ad h{v , K) . {m 

^2E U p X Jyi Jv/k - U(y)/U - y/n \f2E Q U G Xo 



a E = 



2iAU£dk rvc 2 r— — i 2iAU^di 



(Pa(y) + Pe {y)/Z 2 )^/K - U(y)/U - y/Kdy}l = 4 r: 2 f 2 (v,K),(18) 
E*(poXo)t J ^ E*X> 



where y = x/d and Uq is the potential barrier of plane for a straight single crystal and 
the parameter k = U R/(E (3 2 d). This dimensionless parameter one can represent also 
as k = R/Rq, where R = (3 2 E d/U is the characteristic radius of volume reflection [2]. 
Another dimensionless parameter v = E/5E. The parameter v is coupled with the £- 
parameter (see Eq. (9)) by the relation: £ = v — v m = (E — E m )/5E, where E m is the 
nearest local maximum of a transversal energy before the critical point x c {y m = E m /5E). 
Note the functions E(£) and <J E (£) are periodic functions of ^-variable with the period equal 
to 1. 
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The probability e 1 one can presented as sum of two first terms of Taylor series in the 
vicinity of k q : 

£i ( k )=s(k ) + ^(k )k c (--^), (19) 

where k c = UoR c /(Eoj3 2 d) where R c = Eo/£ max is the critical radius of channeling (S max is 
the maximal value of the interplanar electrical field). Thus, k c = Uq/ {£ max d) and k/k c — 
ko/k c = R/R c — Rq/R c (radius Rq corresponds to kq- value). The presentation of a bending 
radius in units of critical radius is convenient for consideration of the process at different 
particle energies. 

From Eq.(9) one can get 

% = -4= T - r - U\ + 7)) ex P (-« + U7(2£ 2 )R, (20) 

where f[, f' 2 is the corresponding derivatives (of fi, f 2 - functions) over /t-parameter. 

Taking into account Eqs.(17 ) and (18) one can represent Eq. (20 ) in the following form: 
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^n(«oK = - t JM, (21) 



where 

JM = ~ f t(t ~ T - J- (- + T)) ex P + UV(2S 2 ))^- (22) 

JO J2 Jl h ?m ^0 J2 

The distribution of scattered particles over the relative transversal energy £' = (E 1 ' — 
E m )/5E follows from integration of Eq.(7): 

, dNg) 1 ride (e'-e-u 2 ^ 

^ ) = ^^ = 72¥io S 6XP 2S^ • (23) 



This function describes the all scattered particles at —00 < £' < 00. The case £' < 
corresponds to particles with £ < i? m and the case £ > 1 corresponds to particles with 
E > E M . It is significant to note that here and below (in similar equations) the values £ m 
and S are functions of the integration variable. Integration of Eq.(23) at condition £ < 
gives Eq.(9) for total number particles, and at condition £' > 1 gives corresponding total 
number of particles with E > Em- 

The distribution of captured particles (due to the additional process, see the trajectory 3 in 
Fig.l) is 

^ '|-^2Hi^!VS^. (25) 



2,71 

Here £ < 1, £ — (E — E m )/5E, £ = (E - E m )/6E, where E and £ are the corre- 
sponding transversal energies (see Fig. 1). Thus, the total normalized on unit number of 
captured particles in this approximation is S\ = S\ a + Eu, where £15 = f^ 00 J I 'bd£" ■ Note our 
consideration of volume capture is valid for small enough bending radii of single crystals. 
The range of applicability of present description we will study below. 
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4 Channeling of the captured particles 



For some time the captured particles move in a planar channel, but due to multiple scatter- 
ing (in the channeling regime) these particles obtain transversal energies which exceed the 
potential barrier and hence the part of them return in underbarrier motion. In this section 
we consider this process on the basis of one-trajectory diffusion approximation. 

For calculations we need to know the losses of transversal energy of a channeled particle. 
For this aim we can use Eqs. (17)-(18). However, we should put the transversal energy 
E m — Ur < E < E m at location of particles in the range x\ < x < x c (Ur = E m — E c h, 
see Fig. 1). At this condition a particle will be in the channeling regime. At first, we find 
the mean and mean square losses of a transversal energy per one period of motion l c (l c is a 
function of a transversal energy). Then we can get approximately that these losses depend 
on the factor l/l c (see Eqs.(10)-(13)), where / is the length of a single crystal. We think at 
l/l c > 1 this approximation is good enough. 

Now we find the distribution function of channeling fraction: 

= 4-/° pyf) + ^(rti ^Ht'-f'-u)'/(gai i( (26) 

V27T J^-U R /8E L, c 

where £ mc = (l/l c )E/SE and S c = (l/l c )a E /5E and the variables £' and £" are the similar 
to analogous ones in Eq(25). Note due to periodicity of the process we wrote £'-value as 
the argument of jF b -function (instead of £ +1). Obviously this equation is valid for small 
enough /-values while the captured particles have transversal energies more than E m — Ur. 

Integration of J- c {£, )-function over £ (from « -U R /5E till 0) gives the normalized 
number of channeling particles £2(0- 

5 Angle distribution of scattered particles 

Relations describing the angle distributions of scattered particles in the process of volume 
reflection were obtained in Ref.[2]. However, these relations do not take into account the 
volume capture. If we want take into account this factor we should divide all particles on 
two sorts. The first sort is the particles which were not captured and the second one is 
the captured particles. The number of particles of the first sort is 1 — S\ and the number 
of particles of the second sort is e±. The summary angle distribution of particles may 
be represented as dNj- / da(a) = dN vr /da(a) + dN vc /da(a), where dN vr /da is the angle 
distribution of pure volume reflected particles (see detail description in Ref. [2]) and dN vc /da 
is the angle distribution of volume captured particles. 

Let us consider the bent single crystal of the finite thickness L. Then the volume captured 
particles we can represent as a sum of the two fraction: the first fraction N c h is the particles 
which enter from a single crystal in the channeling regime (overbarrier motion) and the 
second one A^ e is the particles which transmitted from channeling in underbarrier motion 
(due to dechanneling mechanism). Let us suppose that normalized number particles of the 
first fraction is e 2 , then the normalized number of particles in the second fraction is equal 

to E\ — £ 2 . 

In Ref. [2] the term "physically narrow angle distribution of entering particles" was in- 
troduced. This term means that characteristic size of angle distribution of the entering 
particles" exceeds the angle period 89 = d/{R6) (here 9 is initial angle of particle) and par- 
ticles are uniformly distributed within this period. To this point we consider the multiple 
scattering only in narrow range of transversal coordinate (< d, see Eqs. (17) and (18)). In 



6 



this approximation we can write 



(a) = ^(-<W2 + ^(fe - vMvO- (27) 

here a m is the mean angle of volume reflection, and angle tp — (I — L/2)/R, where tp max — 
L/(2R), and $(x) = 1 or at x > and x < 0), respectively. This equation is valid for 
symmetric case of orientation of a single crystal. Note the case of nonsymmetric orientation 
is considered analogously. 

Fig. lb illustrates the geometry for Eq.(27). Here the broken line COD represents motion 
of a particle due to volume reflection. The lines CE and OD are the directions of particle 
motion before and after process. The angle DOE is the mean angle of volume reflection. AB 
is a tangent line in the critical point. We define the direction along CE-line as zero angle ( 
the initial direction of particle motion). It means that the angle between OD- line (direction 
of volume reflection) and the current direction of channeling flux is approximately equal to 
a = —a m /2 + tp (see Eq.(27)). 

For final result we should take into account multiple scattering in the body of a single 
crystal. In accordance with Ref. [2] we get 

,dN de . r°° . ,dN dei 



("^~( a )) ^ J_ oc p ( ip,ams " , ~d(f ( a ~v)d<P (28) 

where p(tp, a ms ) is the distribution function of multiple scattering with a ms corresponding to 
thickness L/2- (I- L/2)d(l - L/2), I = tpR - L/2 (I - L/2 is the part of particle path in 
the body of single crystal which corresponds to motion in the channeling regime). 

In our model the process of volume capture takes place in the transversal space from xl 
till x3 (see Fig. 1). It corresponds to 5-function distribution of the captured particles. In 
reality, multiple scattering distributed this process in some transversal area and hence in 
some longitudinal space. One can describe the evolution of distribution of volume captured 
particles with the help of the function [2]: 

dN - g i c::p ^- l ) 2 (29) 

where a ca is the mean squared angle of multiple scattering corresponding to approximately 
half of single crystal thickness. 

The density of dechanneling particles as a function of the length / (or the angle tp ) is 
described by the relation as Eqs.(27) and (28) (at the condition a m = 0) but for a ms should 
be taken for a length equal to L/2. It means that we consider particles in the vicinity of 
local variable / and do not taken into account further motion in the body of a single crystal. 
Let us denote this function as dN de0 /dl (or dN de0 /dtp = RdN de0 /dl). 

Now we can describe the behavior of channeling fraction as a function of the length (or 
the angle tp ) by the following equation: 

e 2 (l) = N ch {l) = - dJ ^{z))dz (30) 

Now we can find the angle distribution of channeling fraction. It easy to get the velocity 
(v(t) = x(t)) distribution function (normalized on 1) of channeled particles for one fixed 
transversal energy E: 

^[ () = ,o n 

dv {V> Tc 2 \[3 2 E /R-£(x(t))\ 1 } 
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where r is the half of period of motion for channeled particle and S is the intensity of planar 
electric field. Taking account that 9 v /c we can recalculate this distribution in angle one. 
Knowing the distribution function over the transversal energy of channeled particles (see 
Eq.(26) ) we can find the result distribution of these particles at the exit from a crystal. 



6 Examples of calculations 
6.1 Introduction remarks 

In this section in accordance with obtained equations we present calculations of volume 
reflection process, which taking account the processes of volume capture and channeling 
in bent single crystals. For illustration we selected some energies of proton beams which 
approximately corresponds to energies of well-known accelerators. The main and detailed 
calculations were done for silicon single crystals, but some results were obtained and for 
some another single crystals. For calculation of the potential, electric field, electron and 
atomic densities in silicon we use values of atomic form factors from x-ray experiments. The 
method of calculations one can find in Ref.[8]. Note recent experimental data [7] show that 
the potential from x-ray diffraction gives more precise description of measured characteristics 
of volume reflection for a silicon single crystal than Molier potential. 



6.2 Energy losses over and under barrier particles 

For determination of different distribution functions we should know the mean and mean 
squared energy losses of particles (see Eqs.(17)-(18)). Fig. 2 illustrates these quantities 
for over and underbarrier motion of 400 GeV/c protons. For these calculations we use 
A = HMeV (see Eqs.(10)-(11)). This choice was based on equation for the mean square 
angle of multiple scattering [9] 9 ms = 13.6 [Me V]/ f3cp^/l/X [l + 0.038 In (l/X )\, where p 
is the particle momentum. For a small thickness I I/Xq rs 10~ 3 and more natural law 
6 rns = A/E ^Jl/X the best agreement is at A = HMeV. 

It is easy to recalculate the results of Fig. 2 for any energy of particles (see Eqs.(17)-(18)). 
The universal functions /i and f 2 one can also find from Fig. 2. 



6.3 Simulations of volume capture 

The volume capture probability represents the sum of two terms S\ = S\ a + e^. Initially 
we consider behavior of the first term ei a . Fig. 3 illustrates the calculations of this value 
for different proton energies. One can see that at small values R/R c E\ is approximated 
good enough by linear function. From Eq.(21) we see that value n c dei a /dn is defined by 
J (k ) -integral. Fig. 4 illustrates calculated J-value as a function of k . We see long enough 
plateau beginning from k slightly more than k c . It is convenient to select such k that 
^ia(^i) = 0, where we denote this special value of k - set as k±. It is easy to see that 

= + W (32) 

where kq is any value of k in the range of linearity of S\ a function. According to our 
calculation for the silicon (110) plane K\ 0.13. Now we can wrote the following relation: 

i 

e la = - A ^ Jp - r (- - -) (33) 
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Here we denote the linear approximation of e lo (see Eq.(9)) as E\ a and the plateau value of 

We calculated also value e\b (see Eq.(25) ) for different crystals. According to these 
calculations relation eib/s\ a is approximately equal to 0.41- 0.37 at R/R c = 1.5 — 30, corre- 
spondingly. We think that we can decide this relation approximately constant in this region 
of radii and equal to 0.39. Then we can write the final equation for the probability e\. 

1. 39 AUfJp k ki. . . 

^ ~ 1 " — ) ( 34 ) 

2-^E^S max d2Xi K * «<= 

Analogously, here E\ denotes the linear approximation of e\. 

Fig. 5 illustrates behavior of the relations S\je\ and e la /e la as functions of variable X = 
(k — Ki)/k c . One can see that in the region from X = 0.2 till xs 2.5 these relations are 
constant with a good enough accuracy. Thus, Eq.(34) allows one to describe by universal 
way the probability of volume capture in the area of interest to us. Table 1 contains the 
parameters of some single crystals for calculation of ei-value. 

Fig. 6 illustrates the differential over energy distributions of captured protons which were 
calculated with the help of Eqs. (23) and (25). We see that these distributions become more 
wide with the increasing of the bending radius. 



6.4 Channeling of captured particles 

The channeling of captured particles we illustrate for proton energy 400 GeV and 10 m 
bending radius for the (110) plane of a silicon single crystal. The similar simulations based 
on Monte Carlo calculations for these parameters presented in Ref.[12]. and hence, we can 
compare results which were carried out by different methods. Eq.(26) allow one to calculate 
the flux of particles which were captured in the channeling regime. Fig. 7 shows behavior of 
the flux as a function of the thickness of a crystal. The curves 1 and 2 (see fig. 7b) present 
the speed of dechanneling as a function of variable bending angle for cases without and with 
consideration of multiple scattering in a silicon crystal. In this figure the curve 3 presents 
the distribution function of dechanneling particles on the exit of a single crystal. The total 
number of captured particles is « 0.057. The same number is « 0.064 accordingly to Ref.[12]. 



6.5 Total scattering curve 

Thus, the total scattering curve (with taking account the volume capture process) represents 
the sum of the third curves: 

i) Pure volume reflection curve; 

ii) the curve of the dechanneling fraction; 

iii) the curve of the channeling fraction. 

The resulting curve (for 400 GeV protons and 10 m bending radius of (110) plane silicon 
single crystal) is shown in Figs. 8 and 9 in linear and logarithmic scales, correspondingly. 
The curve of channeling fraction was calculated with accordance of Eq.(31). The semiperiod 
in Eq.(31) was chosen at the condition that in the moment of time t\ the transversal velocity 
is maximal and in the moment of time £2 — ti + t this velocity is minimal. 

One can see that in accordance with the distribution function most particles of channeling 
fraction are groped near critical channeling angle. This fact is due to firstly, the large number 
of particles have high transversal energy and secondly, the time of stay in state with small 
transversal velocity is considerable larger than with the time of state with the large velocity. 
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It should be noted that experimental observation of two peaks of the channeling curve is 
practically impossible because of multiple scattering in detectors. 

The important characteristics of volume reflection of a beam is the efficiency of the 
process which defined as a partition of beam which contained in ±3a vr around of mean 
angle of reflection a m , where a vr is the mean square of angle distribution of a scattered 
beam. The calculated for conditions as in Figs. 8 and 9 the efficiency is close to 0.98. This 
is in agreement with the experimental data [13]. From this result follows that at above 
mentioned definition of the efficiency e the total inefficiency ( 1 — e) should be significantly 
less than the probability of volume capture e\ due to a valuable fraction of particles, which 
quickly return from channeling in underbarier motion (see Fig. 7). 

7 Discussion 

As mentioned above our consideration is valid for a small enough bending radii of single crys- 
tals. It is clear that for large bending radii the energy gap 8E becomes very small. Because 
of this, the particle energy losses can be significantly exceed value 5E, and, consequently, 
the possibility of volume capture will realized for long enough part of trajectory of particles, 
instead of short distance from x\ till x^ (see Fig. 1). 

Let us estimate the area of application of our analytical description. Taking into account 
that distribution of energy losses is described by normal one we can write the following 
condition for validity of obtained in the paper equations: 

<?E,max < 8E/A (35) 

where <JE,max is the maximal value of <je at fixed /t-parameter. For the (110) silicon plane and 
proton energy equal to 400 GeV SE/4[eV] ~ 20/ R[m] and cr E max ~ 1 eV and we got that 
our consideration is valid till R ~ 20 m (or for R/R c ~ 30). The area of validity is extended 
with the increasing of particle energy. It is important to note that the range Rj R c « 10 — 20 
is more preferable for utilization of the volume reflection process on accelerators. Really, at 
these values of radii the mean angle a m of the process is close to maximal one [2] and the 
volume capture process is suppressed in comparing with more large bending radii. 

In this paper for small enough bending radii we got simple relation for probability of 

volume capture S\ in different single crystals. According to this relation S\ is proportional 

_i 

to E 4 (R/R c — Ri/R c ). It means that S\ is a weakly dropping function of particle energy 
at fixed R/R c . For testing of this dependence we carried out Monte Carlo calculations of 
E\ at different radii and energies. The results of these calculations are described by the 
function which is proportional to Eq°' 2 (R/ R c — 0.7). These dependence is close to obtained 
one in this paper. We plan to continue Monte Carlo calculations with the aim of obtaining 
large statistics and hope to find the degree at E which gives the best description of volume 
capture. 

Process of volume capture was investigated also in Ref. [6]. Here the probability of volume 

capture is E\ = R9 c h/L^, where is the dechanneling length, which is proportional to 

_ l. 

particle energy. It means that e\ is proportional to E 2 R/R c . The paper [12] contains the 
criticism of this relation. Authors [12] give the conclusion about incorrectness of the relation 
and try to correct its by introducing of a nuclear dechanneling length. This correction allow 
one to get values of probability close to Monte Carlo calculations (for energy 400 GeV) but 
conserve the previous energy dependence. 
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8 Conclusion 



The main results of our investigation of volume reflection of relativistic particles are: 

1) the simple analytical relations for probability of volume capture are obtained; 

2) propagation of volume captured fraction inside a single crystal is considered; 

3) the analytical method for calculation of the angle scattered distribution of particles, 
taking into account different processes are developed; 

4) presented here equations allow one to find efficiency of the volume reflection process 
but analysys of this problem is beyond of the scope of this paper. 
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Crystal 


Z 


X Q , cm 


d, A 


T T T T 

U ,eV 






«i 


T 

dp 


C {diamond) 


6 


12.14 


1.26 


23.52 


79.70 


0.234 


0.198 


1.71 


Si 


14 


9.38 


1.92 


21.38 


60.0 


0.186 


0.13 


1.49 


Ce 


QO 
6Z 


Z.ZO 


z.UU 


4U.0 ( 


111 AP. 
111.40 


n i sq 


n 1 oq 


i/i/i 

1.44 


W 


74 


0.383 


2.23 


138.34 


490.3 


0.126 


0.0085 


1.29 



Table 1: Parameters of single crystals for calculation of volume capture probability (see 
Eq.(34)). Potentials for diamond and silicon are taken from x-ray measurements [10, 11] and 
for germanium and tungsten Molier potential was used. 
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Figure 1: Scheme of volume reflection and volume capture processes: a) the effec- 
tive potential of bent crystallographic planes, x is the transversal coordinate. Em 
and E m are the transversal energies corresponding to neighbouring local maxima of 
the potential; E c h is one of transversal energies corresponding to local minimum of 
potential; x±,X2,X3 are the transversal coordinates of the local maxima; the curves 
1,2,3 reflect the different possibilities for moving particles. x c is the critical co- 
ordinate for volume reflected particles; b) geometric relations of the process. For 
additional information, see the text. 
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Figure 2: Losses of transversal energies E and cr e (in the (110) silicon plane) for 
volume reflection (a) and for channeling regime (b) as functions of variables (E — 
E m )/5E and (E — E ch )/Uji, correspondingly. The numbers 1 — 5 near curves are 
equal to value of ^-parameter and correspond to ^-quantities. Analogously the 
number V — 5' are equal to value of ^-parameter and correspond to cr^-quantities. 
The values E in (b) are enlarged in three times for clearness. The particle energy is 
equal to 400 GeV. 
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10 20 30 40 50 R/Rc 

Figure 3: £i - value for the (110) silicon plane as a function of relation R/R c for 
different energies 50 GeV (1), 400 GeV (2), 1000 GeV (3), 7000 GeV (4). The 
straight lines near curves are linear approximations (see Eq. (19)). 




Figure 4: J- value as a function of relation R/R c for different conditions. Curves 
1,2,3,4 are for silicon single crystal and proton energies 50, 400, 1000 and 7000 GeV, 
respectively. Curves 5,6,7 are for diamond, germanium and tungsten single crystals 
and proton energy equal to 400 GeV. 
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Figure 5: The relations Si a /e la (a) and e 1 /e 1 (b) as functions of X- value (see 
Number near curves are the same as in Fig. 4. 
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Figure 6: The distributions jF ,jF a ,jF b as functions of ^-parameter for bending radii 
equal to 5, 10, 15 m (curves 1,2,3, respectively). JF = T a + J~b- The curves labeled 
as ' and " correspond to T a -,Tb -distributions and the thick curves to JF one. 



17 




Figure 7: Intensity of channeling fraction (a) and differential particles losses due 
to dechanneling (b) as functions of angle ip. Curves 1 and 2 correspond to con- 
siderations without and with multiple scattering. Curve 3 corresponds to angle 
distribution (over a-angle) dechanneling fraction on the entrance of a single crystal. 
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Figure 8: The resulting distribution of scattered particles at volume reflection taking 
into account influence of volume capture. The energy of protons is equal to 400 GeV, 
bending radius of 2-mm silicon single crystal is equal to 10 m. 
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Figure 9: The curve 1 as the same as in Fig. 8 but y-scale is logarithmic. The curve 
2 is calculated without taking into account the volume capture. 
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